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Abstract 

The one-electron quantum-electrodynamic corrections to the magnetic-dipole transition amplitude be- 
tween the fine-structure levels {ls^2s^2p) '^P^j2 — in boronlike ions are calculated to all orders in aZ. 
The results obtained serve for improving the theoretical accuracy of the lifetime of the {ls'^2s^2p) ^^3/2 
level in boronlike argon. 

PACS numbers: 31.30 Jv, 32.70 Cs 
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I. INTRODUCTION 
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The precision in measurements of decay rates of forbidden transitions has considerably in- 
creased during the last years mmm . The accuracy achieved for the magnetic-dipole 
transition (ls'^2s'^2p) — in B-like Ar became better than one part per thousand 

This experimental precision demands a corresponding increase in the accuracy of the 
theoretical predictions. In a recent work lO] we have calculated the Ml -transition probabilities be- 
tween the fine- structure levels in B- and Be-like ions within the region of nuclear charge numbers 
Z = 16 — 22. In particular, it was found that the theoretical result for the case of Ar^^+ deviates 
from the experimental one by about 3a. 

In Ref. llSl], the relativistic, interelectronic-interaction, and quantum-electrodynamic (QED) cor- 
rections to the Ml -transition amplitude were computed, while experimental values were taken 
for the transition energies. The configuration-interaction method in the Dirac-Fock-Sturm basis 
(CIDFS) was employed in order to evaluate the interelectronic-interaction contribution. Correc- 
tions due to single excitations to the negative-continuum energy states were taken into account in 
the many-electron wave functions. As it is known from ||9|,110|], such corrections may be significant 
in calculations involving operators, which mix large and small components of the wave functions, 
such as the Ml -transition operator. The frequency-dependent term (consult the detailed descrip- 
tion presented in Ref. Illi ') was calculated within perturbation theory to first order in 1/Z. The 
QED correction was obtained within leading order by including the electron anomalous magnetic 
moment (EAMM) in the Ml -transition operator. Uncalculated higher-order QED terms together 
with the experimental errors of the transition energy determine the total uncertainty of the theo- 
retical predictions presented in Ref. |8]. In this work, which is aiming for improvements of the 
evaluation of the radiative effects to the Ml -transition probability in B-like ions, we present the 
exact calculation of the one-electron QED corrections going beyond the EAMM approximation. 

Accordingly, the bound-electron propagator is treated exactly. This approach was already em- 
ployed for the evaluation of the radiative corrections to the decays 2pi/2, 2s, 22^/2 — Is in hydro- 
gen ic ions Il2i and parity nonconserving transitions in neutral Cs and Fr lislfl^ . Besides, in Ref. 
1 15|] the QED corrections to the transition probability between the hyperfine-structure components 
were expressed in terms of the corresponding corrections to the bound-electron g factor. The latter 



ones were calculated to all orders in aZ for the Is and 2s states in Refs. 1116 . 
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Relativistic units (h = c = m = 1) and the Heaviside charge unit [a = e /{Att), e < 0] are 
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used throughout the paper. 



II. BASIC FORMULAS 

The magnetic-dipole transition probability between the one-electron states a and b can be writ- 
ten in the form 

271 

2jVTi ^ ^ ^ 



Wu = -^y,TT\A,,,\\ (1) 



rria mi, M 

where the summation over the photon polarization and the integration over the photon energy 
and angles were carried out. The initial state a is characterized by the angular momentum ja, its 
projection m^, and the energy 6a, while the final state b has the corresponding quantum numbers 
jb, nib, and the energy Sb- The transition amplitude Aim is defined by 



A,M = -\—{b\Tl\a), (2) 



TT 



where uj is the transition energy and denote the spherical components of the Ml -transition 
operator 

T- = ^J.(-)^. (3) 

Here ji denotes the first-order spherical Bessel function and ex. is the Dirac-matrix vector. In 
further calculations we take into account only the first term in the power expansion of ji{ujr), 
since for the case under consideration the transition wavelength is much larger than a typical ion 
size. Accordingly, the Ml-transition operator can be related to the magnetic moment operator 

^ = e [r X ol\/2, 

Ti = ^^[rxa] = ^^/x. (4) 
Utilizing the Wigner-Eckart theorem, the transition probability can be expressed in terms of the 



reduced matrix element of (see, e.g., Ref. 112 ih ). which does not depend on the momentum 



projection M. Therefore, it is sufficient to calculate the transition amplitude for a given projection 
M only. In what follows we take M = and omit the corresponding subscript. 

In this work we focus on the one-loop QED contributions to the transition amplitude beyond the 
EAMM approximation. The vacuum-polarization (VP) and self-energy (SE) corrections, which 
one needs to consider, are diagrammatically depicted in Figs. [T] and |2l respectively. The VP 
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FIG. 1 : Feynman diagrams representing the one-loop vacuum-polarization correction to the transition am- 
plitude. The double line indicates the electron propagating in the external field of the nucleus. The photon 
propagator is represented by the wavy line, while the single photon emission is depicted by the wavy line 
with arrow. 
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FIG. 2: Feynman diagrams representing the one-loop self-energy correction to the transition amplitude. 
Notations are the same as in Fig. Q 

correction corresponding to the diagrams presented in Figs. [Ha) and[TJb) (the electric-loop term) 
has been calculated in the Uehling potential approximation. The diagram depicted in Fig.[Uc) (the 
magnetic-loop term) has the magnetic-interaction insertion into the VP loop. It is known, that the 
contribution of this diagram vanishes in the Uehling approximation. The higher-orders VP terms 
turn out to be rather small and can be neglected. The remaining part of the work is devoted to the 
SE correction. Here we present the formal expressions for the corresponding contributions, which 
were derived at length in Ref. ll22ll . 

The contributions of the diagrams depicted in Figs. Eta) andEtb) are conveniently divided into 
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irreducible and reducible parts. The reducible ("red") contribution of the diagram depicted in 
Fig.|2a) is defined as a part in which the intermediate state energy = €a, respectively €n = £b 
for the diagram presented in Fig.Etb). The irreducible ("irr") part is given by the remainder. The 
latter one can be written in terms of nondiagonal matrix elements of the self-energy operator (see, 
for details, Ref. Im) 



"3^ 



((6|SR(£fc)|5a) + (56|SR(£,)|a)) , 



(5) 



where the perturbations to the wave functions are defined as 



n)(n|/i2|a) 



\n) {n\fi,\b) 



(6) 



n n 

is the renormalized self-energy insertion, which is related to the unrenormalized self-energy 



(7) 



by = E(£) — 7°(5m, where 5m is the mass counterterm. In Eq. dTJ) we use the following 

notations = (1,0;), I{E) = e^a^a'^ D^JE), where D^^{E) is the photon propagator. The 
expression for the reducible part is given by 11221] 

"J 



AA 



red 




{b\^i,\a) {{a\E'{ea)\a) + {b\E'{eb)\b)) , 



(8) 



where S'(ea) = dS{e) / de\e=s^. The contribution of the diagram depicted in Fig.|2tc), known as 
the vertex ("ver") term, is given by the equation ll22ll 

(riil/i^lns) {bn2\I{E)\nia) 



2cu3 i 



ni 712 



(9) 



The irreducible part can be renormalized in the same manner as the ordinary SE correction 



to the energy. This renormalization is well-known and discussed in details in ||23 



2±. The 



ultraviolet divergence in the vertex and reducible contributions can be isolated by expanding the 
bound-electron propagator in terms of the interaction with the field of the nucleus. For our pur- 
poses, it is convenient to decompose the total contribution into zero-, one-, and many -potential 
terms according to the number of interactions with the external field 



(10) 







+ Remainder 



FIG. 3: The potential expansion of tlie vertex diagram. Tiie single line indicates the free-electron propagator 
and the line ending with the cross denotes the interaction with the field of the nucleus. All higher-order 
contributions are contained in the Remainder. 



and 



^^red _ ^^red(O) _|_ ^^red(l) _j_ ^^rcd{2+) 



(11) 



This expansion for the vertex diagram is schematically presented in Fig. |3l In order to remove 
the divergences in the vertex and reducible terms, we consider them together. Combining the 
corresponding parts, we define 



0,1,2+). 



(12) 



It can be shown, that the ultraviolet-divergent terms, which are present in AA^°^^^^ and AA^^'^^^\ 
cancel each other in AA^^^^\ The remaining one- and many -potential terms are ultraviolet finite. 

The zero- and one-potential contributions are evaluated in momentum space, while the many- 
potential term is calculated in the coordinate space employing the partial-wave expansion. The 
scheme for the separate treatment of the one-potential term was also used in previous g factor 
calculations presented in Refs. Ill7 . It improves considerably the convergence of 

the partial-wave expansion in the low and middle Z region. The summarized expressions for the 
reducible correction are similar to those derived for the g factor (see, Ref. |20]). However, for 



the vertex contribution there are some principal differences. In contrast to Eq. (4) of Ref. 
two different energies £a and enter the denominator of formula ^ and the matrix element 
{bn2\I{E)\nia) depends on both states a and b. Taking these differences into account, we derive 
the corresponding formulas for the AA"^"^^^^ term in the Appendix. The derivation of the formu- 
las for the one-potential vertex contribution is somewhat more complicated. However, taking the 
energy to be the same for both electron propagators (e.g., Ea), the expressions for AA^°^^^^ can be 
obtained in the same manner as for the g factor [12011 . The remaining many-potential term can be 
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evaluated by the point-by-point subtraction of the corresponding zero- and one-potential contribu- 
tions in the coordinate space. Consistently, we subtract the one-potential vertex contribution with 
the same energy variable in the electron propagators as it is taken in the AA^^^^^^ term calculated 
in the momentum space. Furthermore, the term with = €b and = Sa in Eq. ® has an 
infrared divergence, which is canceled by the corresponding term of the reducible contribution. 

For the numerical evaluation we employ the finite-basis-set method for the Dirac equation con- 
structed via the dual kinetic balance approach 11261] . The summation of the partial-wave expansion 
was performed up to l^maxl = 10, while the remaining tail was approximated by a least- square 
inverse-polynomial fitting. 



III. RESULTS AND DISCUSSION 



The one-loop QED corrections beyond the electron anomalous magnetic moment approxima- 
tion are conveniently expressed in terms of the correction S, which is defined through 

a 

AAqeD = Anr i2Ke + -6) , (13) 

where A^^ is the nonrelativistic transition amplitude and 

Ke = 0.328 478 965... (-) +••• (14) 

27r \TT/ 

represents the electron anomalous magnetic factor corresponding to the EAMM term. In Table 
I we present our results for the one-electron SE correction. The VP term calculated within the 
Uehling approximation has been found to be negligible. The various contributions corresponding 
to the SE corrections to the transition amplitude are given in Table I. The one- and many-potential 
terms are represented as the sum = S^^W + As one can see from this table, the oc- 

curring cancellation reduces the total value for the correction 5 by an order of magnitude compared 
to the individual terms. Most serious computational difficulties arise from the extrapolation of the 
partial-wave expansion of the many-potential term. In order to estimate the error, we perform a 
second evaluation of without separating out the one-potential term. The difference between 

the results of both calculations is taken for the uncertainty. 

The obtained results allow for an improvement of the theoretical values of the Ml -transition 
probabilities between the states (ls^2s^2j9) ^^3/2 — in B-like ions presented in |l8|]. In Ta- 

ble II the results of the CIDFS calculation without the QED term W'^ the improved radiative 
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TABLE I: Individual contributions to the one-electron self-energy correction expressed in terms of the 
various 6 corrections, defined by equation ( fOl . Numbers in parenthesis represent error in the last digit. 



z 


^irr 


^vr(O) 


jvr(l+) 


6 


16 


0.0177 


-0.0124 


-0.0065(1) 


-0.0012(1) 


17 


0.0195 


-0.0136 


-0.0075(1) 


-0.0016(1) 


18 


0.0213 


-0.0148 


-0.0084(1) 


-0.0019(1) 


19 


0.0232 


-0.0160 


-0.0094(1) 


-0.0022(1) 


20 


0.0252 


-0.0172 


-0.0105(1) 


-0.0025(1) 


21 


0.0272 


-0.0185 


-0.0116(1) 


-0.0029(1) 


22 


0.0292 


-0.0197 


-0.0128(1) 


-0.0033(1) 



correction AWqed, the total values of the transition probability Wtotah and the lifetimeTiotai are 
compiled. For the transition energies we used the experimental values from Refs. IitHisI . They 
are presented in the second column of Table II. For S^^+, Cl^^+, K^^+, and Ti^^+ the uncertainties 
of the values of PVtotai and rtotai are determined by the errors in the experimental transition ener- 
gies. Thus the accuracy of the total results for these ions has not been improved in comparison 
with the values presented in Ref. [i8|fl. However, for Ar^^+ the experimental value for the transition 
energy is known with high precision l28i and the uncertainty of the predictions obtained in Ref. ^ 
is determined by the uncalculated QED terms beyond the EAMM approximation. The present cal- 
culation of the radiative correction AVFqed has improved the accuracy of the total values of the 
transition probability VFtotai and the lifetime rtotai for Ar^'^^ by an order of magnitude. Further- 
more, for argon ion we have added the probability of the electric-quadrupole mode We2 = 0.00194 
s~^ calculated in Coulomb gauge in Ref. [29] to the total decay rate and lifetime. 

In Table II, we also compare our total results with correspondingexperimental data. The dis- 
agreement with the most accurate experimental value for Ar^'^+ Ba, LZD can be stated. The reason 
for this discrepancy is still unclear for us. 
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TABLE II: The decay rates W [s^^] of the transition (15^25^2^) ^Pg/a - and the lifetime r [ms] 

of the (ls^2s^2p) state in B-Uke ions. The transition energies are given in cm^^. The values of 
are taken from Ref. ISfl. Results of the present work AWqed> ^totai> and Ttotai are given in columns 4 to 
6. For comparison, the experimental values Texpt are presented in the last column. Numbers in parenthesis 
denote the estimated uncertainty. 



Ions 


Energy and Ref. 


[8J 


AWqed 


W^total 


"n^otal 


Tcxpt and Ref. 


511+ 


13135(1) [27] 


20.34481 


0.09444 


20.439(5) 


48.93(1) 






17408(20) [27] 


47.34975 


0.21972 


47.57(16) 


21.02(7) 


21.2(6) [5] 














21.1(5) [5J 


Ari3+ 


22656.22(1) [28J 


104.36006 


0.48419 


104.846(3) 


9.5378(3) 


9.12(18) [I] 














9.70(15) \J} 














9.573(4)(5) r6J 


Ki4+ 


29006(25) \7T} 


218.9394 


1.0156 


220.0(6) 


4.546(12) 


4.47(10) \£ 


Til7+ 


56243(4) [27] 


1594.714 


7.393 


1602.1(5) 


0.6242(2) 


0.627(10) \2} 



02-17574) and by INTAS-GSI (Grant No. 03-54-3604). A.V.V. and G.R acknowledge financial 
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Appendix: Zero-potential vertex term 

Let us start from the momentum representation of the transition amplitude Aim- Referring to 
Eqs. dSJ) and Q it can be written as 

AiM = -^'^teJ -|-^^,(p) [a X Vp,(5=^(p-p')]^,^a(p'), (15) 

where the gradient Vp/ acts only on the 6 function. In order to obtain the zero-potential ver- 
tex term AA^°^^^\ we substitute ct by the renormalized part of the free-electron vertex operator 

^^ver(o) ^ -^^e I ^^Mp) [rR(e,,p,5,,p') x Vp,5^(p - p')] , ^a(p') • (16) 
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In contrast to the g factor (see Eq. (15) of Ref. oKjH ). the wave functions of the initial (a) and final 
(6) states enter into Eq. (fT6b and rR(eb, p, Ea, p') has different energy arguments. Integrating by 
parts and performing the integration over p' yields 



^^ver(O) 




dp 

{2ny 



dp 

{2ny 



(17) 



where 



{eb,ea,p) = [Vp/ X rR(£b,p,e„,p')]Jp'=p 



(18) 



The right side of Eq. dTTl) is naturally divided into two parts Ayl™'(°)'^ and AA'"''^")'^. Starting 
with the first one, it is convenient to represent the function E(eb, £a, p) in the form 

\ A ■ [ d^/i; 1 i)- /k + m . _ , jk + m „ 

p) = taay X Vp|. ^^,_^^,_ ^^, 7 (19) 

with p = {sb, p), p' = {ea, p), and /jj = p^7^, respectively. Using the commutation identity for 
the 7 matrices, we get 

„ a C d^k 1 

^(£6,£a,P) 



Am'^ J k'^ [{p — ky — m?\ [{p' — k^ — m^] 



[7X (p-k)],7-}. (20) 



{p' — ky — 

Expressing the integration over the loop momenta k in terms of the integrals over the Feynman 
parameters, one can derive the formula 



a 



12) 



+2p\Au - A21) + 2p'\Ai2 - A22) - Ap%A23\ [7 X p].} , (21) 

which coincides with the corresponding equation in calculations of the g factor, if one considers 
Ea = £b- Here the Feynman integrals are determined as 



Co 



{yp+{^-y)p'f 



In X) 



(22) 



C12 



{yp+{l- y)p') 



1-y 



;i-Fln X), 



(23) 
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An 







-f 

Jo 



x{l - x){l - y) 



dxdy 



An 

Al2 



yi_y 



(24) 



(25) 



(26) 



^21 



122 



^^23 y 



^0 



a;^(l — a;)(l — y) 

Z2 



y 

(1-2/)' 



(27) 



and 



1 - yp' - (1 - 2/)p'' 



{yp+{l-y)p'y 
Z = + (1 - + 1 - yp2 _ (1 _ ^)p/2_ 

To carry out the angular integration for the transition under consideration 2p3/2 
employ the following results for basic integrals (/i = 1/2) 

-2V2/3 for Ki = 1, K2 = -2 , 
for = —1, K2 — 2 , 



(28) 
(29) 

2pi/2 we 
(30) 



dJ^pX«iu(p)[o- X p] 



(31) 



—V2/3i for = —1, = — 2 , 
-\/2/3i for = 1, ^2 = 2 , 

where p = p/|p|, X^fiip) is the spherical spinor, and cr denotes the vector of Pauli matrices 
Finally, for the first part AA^^^^'^^'^ we obtain 



+Pr 



[eaEb - Pr)iAo -Ai- 4A23) - {Ao + 3Ai) + 2{el-pl 



x(An - A21) + 2(£2 - - A22)J {gbfa - fb9a) 

-plisa - eb){Ao - Ai){gbga - hfa) \ , 



(32) 
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where j9r = \p\, QaiPr) and fa{Pr) are the upper and lower radial components of the wave function 
in the momentum representation, respectively. 



The second term A/ 



ver(0),2 



can be calculated similarly. Using the expression for the free- 
electron vertex function ll25h and employing in addition to Eq. (BTT ) the following angular integrals 

r , f V2/3i for Ki = 1, K2 = -2, 

dn^xl,mP^^nM^'Ut>) = \ (33) 
J t for = -1, K2 = 2, 

f -2^/2/3 i for Ki = -1, K2 = -2 , 

d%xl,m^y<VnM^.M = \ ' ^^^^ 

y \/2% for = 1, = 2 , 
where Vf^p is the angular part of the gradient, we have 

3 2 3 2 

+ —9bfa fb9a) - PrD{£a " Eb){fbfa ' QbOa + —fbfa gbQa) 

Pr Pr Pr Pr 

+ {e,B + 2e,D + e^C + AD)g,ga+Pr{B + 2D + C)hg}^ , (35) 

where g'^iPr) = dga{pr)/dpr, fa{Pr) = dfa{Pr)/dpr- Here the set of coefficients are expressed in 
terms of the Feynman integrals as 

A = C2A-2 + {el- pl)Cn + {el - pI)Cu + ^SaSb - Pr){Co + Cn + Cu) 

-2Co + Cn + Cu , (36) 

B = -4(Co + 2Cn + Cu + C21 + C23) , (37) 

C = -4(Co + Cn + 2Ci2 + C22 + C23) , (38) 

D = 2(Co + Cn + Ci2), (39) 



and 



C24 = - / dy \n{y'^{ea - - y{ea - e^f + 1) 
Jo 



(40) 



The total result for the zero-potential vertex contribution is the sum of the corresponding terms 
from Eqs. (l32ll and (l35l). 
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